Charge-order- induced sharp Raman peak in Sri 4 Cu 2 404i 
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In the two-leg 5=1/2 ladders of Sri4Cu24 04i a modulation of the exchange coupling arises from the 
charge order within the other structural element, the CuC>2 chains. In general, breaking translational 
invariance by modulation causes gaps within the dispersion of elementary excitations. We show that 
the gap induced by the charge order can drastically change the magnetic Raman spectrum leading 
to the sharp peak observed in Sri4Cu24 04i. This sharp Raman line gives insight in the charge-order 
periodicity and hence in the distribution of carriers. The much broader spectrum of LaeCa8Cu24 04i 
reflects the response of an undoped ladder in the absence of charge order. 

PACS numbers: 75.40. Gb, 75.10.Jm, 74. 25. Ha, 75.50.Ee 
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Understanding the complex interplay of spin and 
charge degrees of freedom in doped quantum-spin sys- 
tems is a key issue in condensed matter physics. 
This interplay governs in particular the physics of 
the planar high-T c superconducting cuprates. In 
the telephone-number compounds A14CU24O41 it gives 
rise to a variety of interesting ground states. In 
La6CagCu2404i the ladders form an insulating spin liq- 
uid [1]. Sri4_ a; Ca a; Cu2404i becomes superconducting 
under external pressure for x ^ 11.5 [2,3] whereas an in- 
sulating charge-ordered state is favored for x ;$5 [1,4-9], 
although the average copper valence does not depend on 
x. The different properties are usually attributed to the 
different distribution of charges between the two subsys- 
tems [10,11], CU2O3 two-leg ladders and Cu02 chains 
[12]. But the interesting physics linked to the differing 
periodicity of ladders and chains [13] is usually neglected. 

Raman scattering offers a powerful tool to examine the 
spectral density of magnetic excitations and thus pro- 
vides important information on the kinetics and on the 
interactions of the elementary excitations. In the un- 
doped spin liquid La6CagCu2404i, the ladders show a 
very broad two-triplet Raman line with slightly differ- 
ent peak positions for leg-leg and rung-rung polarization 
[14]. This agrees very well with theoretical results [15]. 
A very different line shape, however, is found in charge- 
ordered Sri 4 Cu240 4 i, which shows a peculiar sharp peak 
that is observed at the same frequency in both polariza- 
tions [14,16]. This sharp response poses a challenge to 
the understanding of the cuprate ladders and offers the 
opportunity to study the interplay of spin and charge 
degrees of freedom in this fascinating system. 

The sharpness of the Raman peak in Sri4Cu2404i is 
in strong contrast to the very broad two-magnon Raman 
line observed in the undoped two-dimensional (2D) high- 
Tc cuprates, which is still the subject of controversial dis- 
cussions. Gozar et al. [16] argued that the observation of 
a very sharp two-triplet Raman line in a one-dimensional 



(ID) spin liquid suggests that the large width found in 
2D cannot be attributed to quantum fluctuations. 

Here, we challenge this point of view by providing a 
clear explanation for the Raman data of Sri4Cu2404i. 
We propose that the charge-order superstructure gives 
rise to a modulation of the exchange coupling along the 
ladders. The concomitant backfolding of the dispersion of 
the elementary triplet opens gaps at the crossing points. 
This in turn can have a drastic effect on the Raman line 
shape, which we calculate using continuous unitary trans- 
formations (CUTs) [17,15]. The high resolution accessi- 
ble by the CUT approach is decisive to account for the 
very narrow peak we are aiming at. Our results with 
and without charge order excellently describe the Raman 
data of Sr!4Cu2404i and LagCagCu2404i, respectively. 

For zero hole doping, the minimal model for the mag- 
netic properties of the S — 1/2 two-leg ladders in 
A14CU24O41 is an antiferromagnetic Heisenberg Hamil- 
tonian plus a cyclic four-spin exchange term H cyc [18-20] 

H = Jj_ ^ Si,iS^2 + J|| ^ S, lT Si+i !T + H cyc (la) 

H cyc = Jcyc ^(i,l),(i,2),(i+l,2),(i+l,l) ( lb ) 



-^1234 = 

(S 1 S 2 )(S 3 S 4 ) + (S 1 S4)(S 2 S3) 



(lc) 



(S 1 S 3 )(S 2 S 4 



where i denotes the rungs and r 6 {1,2} the legs. The 
exchange couplings along the rungs and along the legs are 
denoted by J± and Jm , resp.. There is also another way to 
include the leading four-spin exchange term by cyclic per- 
mutations [21,20] which differs in certain two-spin terms 
from Eq. (1) [21]. Both Hamiltonians are identical except 
for couplings along the diagonals if Jj_ and J|| are suit- 
ably redefined [22]. Here, we use Hamiltonian (1) since 
it is established that the four-spin terms are the signifi- 
cant ones [24]. The exchange parameters determined in 
Ref. [20] for La5.2Ca8.sCu2404i correspond in our nota- 
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FIG. 1. Scheme [23] of the superstructure along the chains 
and the ladders (c axis). 10 chain units (top row) match 7 
ladder units (bottom row) inducing a modulation in the lad- 
ders with wave vector Qs = ciadder/cchain = 10/7 = 3/7 + 1 
(in r.l.u. of the ladder) [13]. In Sri4Cu2404i the charge order 
(CO) implies an additional superstructure with Qco = 2/10 
(in r.l.u. of the chain) [4-7], corresponding to a peri- 
odicity of 5 • c c hain. It is visualized (middle row) as 
two units of "spin-hole-spin-hole-hole" per 7 rungs (grey 
squares denote the six holes per f.u.). This superstruc- 
ture induces an additional modulation in the ladder with 
Qco = ciaddor/ (5 -Cchain) = 2/7 (in r.l.u. of the ladder). 



tion [22] to J||/Jl = 1.22 ± 0.05, J cyc /Jj_ = 0.21 ± 0.03 
and J± = 1150 ± 100 cm -1 . Using the average values, 
we find within 5% the same values for the spin gap and 
for the two-triplet bound-state energies as in Ref. [20]. 

The exchange coupling in the CuC>2 chains is much 
weaker than in the ladders since it is mediated via Cu- 
O-Cu bonds with an angle close to 90° . Thus the chains 
do not contribute directly to the Raman line at w 3000 
cm -1 . But the presence of the chains gives rise to 7 
inequivalent ladder rungs per formula unit (f.u.) and 
thereby induces a modulation in the ladders (see Fig. 
1). This modulation is characterized by the wave vector 
Qs = 10/7 = 3/7+1 in reciprocal lattice units (r.l.u.) 
of the ladder. In the magnetic subsystem of the spins on 
the Cu sites of the ladder, wave vectors are only mean- 
ingful modulo unity so that Qs = 10/7 and Qs = 3/7 
are equivalent. The additional modulation induced by 
the charge ordering on the chains below Xco ~ 200 K 
has the wave vector Qco = 2/7 [4-7], see Fig. 1. 

Now we estimate the amplitude of the exchange mod- 
ulation with Qs- The Cu-0 distances within the lad- 
der are hardly affected by the modulation; the main ef- 
fect is a shift of the O ions perpendicular to the Cu-O- 
Cu bonds [13]. Hence the electronic hopping elements 
t p d can safely be considered constant. The exchange 
couplings are modified by the induced variation of the 
charge-transfer energy A ct , i.e. the variation of the en- 
ergy difference between holes on Cu and on O. We com- 
puted the variation of A ct in a point-charge model with 
stoichiometric valencies except for the chain oxygen with 
q = — 1.7e to account for the holes [10,11]. The calcu- 
lation uses Ewald sums so that the results pertain to 
the infinite system. The relative changes A Jj J are esti- 
mated in leading order [24]. Assuming structurally un- 



modulated chains and ladders, we find a negligible effect 
of the chains on the exchange couplings of the ladder of 
| A J I J\ £ 10~ 6 . However, the modulated positions at 300 
K [13] yield 



J| M = J|| [1 + 0.05 cos(27rf (i + i))] 



(2a) 



J±,i = J±[l - 0.10 sin(27rfi) +0.05 cos(27r|(i+3))] (2b) 

with phase accuracy \Ai\ J$ 0.1 [23] where i counts the 
leg- or the rung-bonds. The term with 2Qs = 6/7 de- 
notes the second harmonic; overtones with amplitude J$ 
1% are omitted. The amplitudes in Eq. (2) show that the 
induced modulation of the couplings is indeed sizeable. 

We expect that the effects of the charge order occurring 
below Too ~ 200 K are of similar size. Without detailed 
information on the structure at T <C Tco, only an esti- 
mate is possible. We assume a charge modulation on the 
chain oxygen of Aq(j) = — 0.2e cos(27r^(j + |)), where j 
counts the chain O sites, and the periodicity 5c c h a in and 
the phase are established experimentally [4-7] (cf. middle 
row of Fig. 1). This yields an additional modulation 



AJ !M = 0.16J||COs(27r|i) 



(3) 



with phase accuracy of |Ai| ;$ 0.1 [23]. Thus the modu- 
lation induced by the charge order with Qco = 2/7 (in 
r.l.u. of the ladder) is indeed significant. 

Now, we investigate the effects of modulations on mag- 
netic Raman scattering which probes the excitations with 
zero momentum and zero spin. At T—0 the Raman re- 
sponse I(w) is given by the retarded resolvent 



I(uj) = ^- lim lm(0 



i-H+Eo+iS 



R 



(4) 



The observables R lun s (_R lc s) for magnetic light scatter- 
ing in rung-rung (leg-leg) polarization are given in Ref. 
[15]. We focus on the dominant two-triplet contribution. 
A CUT is employed to map the Hamiltonian H to an 
effective Hamiltonian H e g which conserves the number 
of rung-triplets [17,25,15]. The ground state of H e s is 
the rung-triplet vacuum. The observable R in I(to) is 
mapped by the same unitary transformation to an effec- 
tive observable i? e ff ■ 

The CUT is implemented perturbatively in J|| / J± . We 
compute H e g and R c ff to order n > 10. A calcula- 
tion in order n accounts for hopping and interaction pro- 
cesses extending over a distance of n rungs. The result- 
ing plain series are represented in terms of the variable 
1 — A/ ( J|| + J±) [26]. Then standard Pade approximants 
[27] yield reliable results up to J\\/J± ~ 1 — 1.5 depend- 
ing on the value of J cyc /^_L- Consistency checks were 
carried out by extrapolating the involved quantities be- 
fore and after Fourier transforms. In case of inconclusive 
approximants the bare truncated series are used. We es- 
timate the overall accuracy to be «5%. The Raman line 
shape is finally calculated as continued fraction by tridi- 
agonalization of the effective two-triplet Hamiltonian in a 
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FIG. 2. Dispersion of the elementary triplets with (thick 
line) and without (thin line) a modulation of «2/7=0.15 (cf. 
Eq. 5) with Qco — 2/7 (arrows). Some higher-order contri- 
butions are denoted by dashed arrows. 

mixed representation using the total momentum and the 
real-space distance. So the total momentum is sharply 
defined. No finite-size effects appear. This ensures a 
particularly high resolution in momentum and in energy 
necessary to account for a very sharp feature. 

The modulation is included on the level of the effective 
model, i.e. after the CUT. This is no serious caveat since 
a microscopic calculation is not available. The leading- 
order effect of J|| is to enable the elementary triplets to 
hop from rung to rung by a nearest-neighbor hopping 
element t\ oc Jm and to induce a nearest-neighbor in- 
teraction Wi oc Jm . So the most straightforward way to 
account for the modulation of J|| as given in Eqs. (2a) 
and (3) is to modulate t x and W\, 

ti oc w\ oc J|| ■ [ 1 + ctQ cos(2irQi) } . (5) 

Q=2/7, 3/7, 6/7 

Since we focus on the effect of the charge order (Eq. 3), 
the modulation of J± as given in Eq. (2b) is neglected. 

We use the parameters fixed for La5.2Ca8.sCu2404i 
in Refs. [20,22], J||/Jj_=1.25 and J cyc /Ji=0.2. Fig. 2 
shows the dispersion with and without a 15% modula- 
tion with wave vector QcOj i- e - 0:2/7=0. 15 in Eq. (5). 
Clearly, sizeable gaps open wherever Qco links equal 
energies uj(k) = ui(k + Qco) of the unmodulated lad- 
der. Smaller gaps open for higher-order processes, e.g. 
for uj(k) = uj(k + 2Qco). Thus the energies at which 
gaps open depend decisively on the wave vector of the 
modulation. 

The Raman response of an undoped and unmodulated 
ladder is very broad (bottom panel of Fig. 3 and Ref. 
[15]), in good agreement with data of La6CasCu2404i 
[14] (middle panel of Fig. 3). The excellent description 
of the peak position obtained for the parameter set of Ref. 
[20] given above corroborates these parameters. What is 
the effect of a modulation on the Raman line shape? The 
occurrence of gaps implies prominent peaks (van Hove 
singularities) in the density of states (DOS) and hence in 
the Raman line shape. Since Raman scattering measures 
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FIG. 3. Thick lines: Raman line shapes in leg-leg polar- 
ization for J\\/J± = 1.25, J cy c/J± =0.2 and J± = 1100 cm -1 
[20,22] without modulation (bottom), with the structural 
modulation Qs—S/7 and 6/7 appropriate for La6CagCu2404i 
(middle) and with the additional charge-order modulation 
of Sri4Cu2404i (top). Thin lines: Raman data from Ref. 
[14], T=20K. A modulation-induced gap in the dispersion 
at ujg (see Fig. 2) causes a Raman feature at 2u) g . Addi- 
tional features arise due to backfolding, e.g. the small peak at 
2200 cm -1 corresponds to the S=0 two-triplet bound state at 
k=2/7 [25]. 

excitations with total momentum fc tot =0, the two-triplet 
response reflects the excitation of two triplets with mo- 
menta k 2 — —ki and energies oj(ki) — uj(k 2 ). A gap at 
uj g thus causes a corresponding feature in the Raman 
line at 2u g . For the structural wave vectors Qs—'i/l 
and 6/7 these effects are rather small (middle panel of 
Fig. 3). But a drastic change of the line shape appears 
if (and only if) 2uo g coincides with the broad peak of 
the unmodulated ladder, since then the opening of the 
gap implies a redistribution of a large part of the spec- 
tral weight. For the relevant exchange couplings we find 
that 2 • w(7r/ci a ddcr ± Qco/2) « 2.8 J±_ w 3100 cm" 1 is 
slightly above the Raman peak of the unmodulated lad- 
der. Hence the charge-order modulation piles up a large 
part of the high-frequency weight on top of the peak, 
giving rise to a very sharp feature which agrees very well 
with the data of Sri4Cu2404i [14,16] (top panel of Fig. 
3). Since both the exchange constants in Sri4Cu2404i 
and the calculations are accurate within a few percent, 
only a semi-quantitative analysis is possible which shows 
the principal mechanism. The remaining uncertainties 
may imply that also a smaller value of a 2 /j is sufficient 
to produce the sharp feature. 

The good agreement between the experimental data 
and the theoretical result, based on the independently 
determined couplings [20,22] and the wave vector of the 
charge order, corroborates our interpretation. Another 
argument stems from the polarization dependence. For 
Jcyc > 0, the peak positions for leg-leg and rung-rung 
polarization should be different [15]. This is indeed 
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the case in LaeCasCu^Cn [14], but the sharp peak in 
Sri4Cu2404i is found at w 3000 cm -1 in both polariza- 
tions [14,16]. In the scenario of the modulation-induced 
gaps the peak position is determined by the position and 
the size of the gap, since the peak is primarily a DOS 
effect. Hence the coincidence of the peak positions in 
both polarizations despite J cyc > supports our scenario. 
Moreover, also the spectra of Sri4Cu2404i at elevated 
temperatures are explained. The charge order melts at 
Tco ~ 200 K. Indeed, the very sharp Raman line is ob- 
served only below Tco [14]. For T > Tco, the peak po- 
sitions are different for the two polarizations [14], which 
is expected for J cyc w 0.2J± at a 2 /7 = 0. 

An alternative explanation of the sharp peak in terms 
of bound states is unlikely. There is no bound state 
within the broad Raman continuum of the undistorted, 
undoped ladder [15]. But how about finite doping? At 
300 K, about 90 % of the doped carriers reside in the 
chains [10]. At low temperatures the distribution of 
holes is not yet settled experimentally. Theoretically, the 
Madelung potentials indicate that all the holes reside in 
the chains [11]. This is corroborated by the observation 
of the periodicity 5c c h a in [4-7] for the charge order in the 
chains. In a ID fermionic system it is natural to view the 
charge order as an effect of the 2/cp instability. So we are 
led to conclude that 2fcp = 2/10 (in r.l.u. of the chain), 
which implies that there are nj, + n-f = 4/10 electrons 
per site or 6 holes per f.u.. This further supports our 
assumption that at T sa all holes reside on the chains. 

A finite hole concentration on the ladders cannot be 
ruled out completely These charges would be pinned 
in a commensurate charge-density wave (CDW) at low 
temperatures by the electrostatic potential of the charge 
order on the chains. Clearly, such a CDW would also 
induce strong modulations. But it remains unclear where 
the peculiar periodicity stems from if 2kp ^ 2/10. 

A small amount of impurity holes cannot explain the 
sharp Raman peak. Below Tco, Sri4Cu2404i is insulat- 
ing, i.e. all charge carriers are localized. The local charge 
degrees of freedom may couple to the magnetic ones, but 
due to the local character the whole Brillouin zone would 
be involved implying a broad energy distribution, at odds 
with experiment [14,16]. 

Our results clearly call for several experimental ver- 
ifications. Neutron-scattering experiments could clarify 
the presence and the size of gaps in the dispersion. Low- 
temperature investigations of the structure would help to 
improve our understanding of the charge-ordered state. 
Low-temperature x-ray absorption data are required to 
determine the hole density in the ladders. Raman stud- 
ies as a function of Ca concentration and temperature 
could verify that the features explained here are indeed 
connected to the occurrence of the charge-ordered state. 
Then, the peak position offers a sensitive tool to deter- 
mine the modulation wave vector Qco- 

In conclusion, the modulation of the exchange coupling 



in the charge-ordered state of Sri4Cu2404i can explain 
the peculiar Raman data. The induced gap redistributes 
a large part of the spectral weight, giving rise to a sharp 
Raman peak. A comparison with the 2D cuprates is not 
appropriate. Strong quantum fluctuations are still the 
most likely candidate to explain their very broad Raman 
line. 
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